Draft version August 6, 2010 

Preprint typeset using I^TgX style emulateapj v. 2/16/10 



O 
(N 

< 

o 
u 

6 

> 

in 

vn 

o 
o 



X 



REVERSAL OF FORTUNE: CONFIRMATION OF AN 
INCREASING STAR FORMATION-DENSITY RELATION IN A CLUSTER AT Z = L62i 

Kim-Vy H. Tran^-\ Casey Papovich^, Amelie Saintonge'', Mark Brodwin^, James S. Dunlop'', Duncan Farrah', Keely 
D. Finkelstein', Steven L. Finkelstein-, Jennifer Lotz*, Ross J. McLure*, IvelinaMomcheva', & Christopher N. A. 

Willmer'" 

Draft version August 6, 2010 

ABSTRACT 

We measure the rest-frame colors (dust-corrected), infrared luminosities, star formation rates, and stellar 
masses of 92 galaxies in a Spitzer-selected cluster at z = 1 .62. By fitting spectral energy distributions (SEDs) to 
10-band photometry (0.4/im< X„bs < 8/im) and measuring 24/im fluxes for the 12 spectroscopically confirmed 
and 80 photometrically selected members, we discover an exceptionally high level of star formation in the 
cluster core of ^ 1700 M0 yr"' Mpc"^. The cluster galaxies define a strong blue sequence in (U - V) color and 
span a range in color We identify 17 members with L/r > 10" Lq, and these IR luminous members follow 
the same trend of increasing star formation with stellar mass that is observed in the field at z ^ 2. Using rates 
derived from both the 24/im imaging and SED fitting, we find that the relative fraction of star-forming members 
triples from the lowest to highest galaxy density regions, e.g. the IR luminous fraction increases from ^ 8% at 
E ~ 10 gal Mpc"^ to ^ 25% at S > 100 gal Mpc"^. The observed increase is a reversal of the well-documented 
trend at z < 1 and signals that we have reached the epoch when massive cluster galaxies are still forming a 
substantial fraction of their stars. 



Subject headings: galaxies: evolution - galaxies: starburst - 
-infrared: galaxies 



galaxies: clusters: individual (CIG J02 18. 3 -05 10) 



1. INTRODUCTION 

A well-established observational hallmark of how galax- 
ies evolve as a function of environment is the star formation- 
galaxy density relation. A plethora of studies utilizing multi- 
wavelength tracers of activity have shown that star forma- 
tion universally decreases with increasing galaxy density 
at z < 1 (e .g. Hashimoto etal.'l99i lEUingson et al] I2OO1I 
iGomez eFa l. 2003; Patel et al. 2009i). In particular, the cores 
of massive galaxy clusters are galaxy graveyards full of mas- 
sive spheroidal systems that are dominated by old stellar pop- 
ulations. However, as we approach the epoch when these 
quiescent behemoths sh ould be forming t he bulk of their 
stars (z > 2; IvTn Dokku m'etan iT998: J0rge nsenet all 120061: 
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iRettura et al.l201C)l) , the star formation-density relation should 
weaken and possibly reverse. Identifying when star forma- 
tion is quenched as a function of galaxy mass and envi- 
ronment provides strong constraints on galaxy models (e.g. 
Kauffmann et al. 1993; Hopkins et al. 2008), i.e. is individ- 
ual galaxy mass or the larger scale environment the primary 
driver of evolution? 

Observations of field galaxies at z '^ 1 indicate that the star 
formation-density relation turns over at this epoch such that 
there is an enhance ment of activity in the highest densit y 
regions of the field (Elb az et al.1 120071: ICooperet all l2008j. 
Studies also find an excess of dust-obscure d star forma tion in 
group environments at z < 1 JKovama et al. 2008; Tra n et al.l 
I2009t iGallazzi et al.l 120091) , and very recent results suggest 
that star formation may be enhanced i n the significantly r icher 
core of a galaxy cluster at z ^ 146 (iHilton et al.ll201(> '). As 
cluster surveys push to higher redshifts (z > 1) and thus closer 
to the epoch when massive galaxies are forming their stars, 
variations in age will become evident in, e.g. a larger scatter 
in color, and robust star formation rates should reveal increas- 
ing levels of activity even in the cluster cores. 

We report here the first confirmation of increasing star for- 
mation activity with increasing galaxy density observed in 
cluster CIG J021 8. 3-05101 , a Spi tzer-selected galaxy cluster 
at z = 1.62 (Pa povich et al.ll20Tol hereafter PaplO). We use 
cosmological parameters fi,,, = 0.3, A = 0.7, and H = 70 km 
s"' Mpc"' throughout the paper; at z = 1.62, this corresponds 
to an angular scale 0.5 Mpc arcmin"'. All magnitudes are in 
the AB system. 

2. MULTI- WAVELENGTH DATA 

CIG J0218. 3-0510 has a wealth of multi-wavelength imag- 
ing data that in cludes BRi'z' imagin g from the Subaru-XMM 
Deep Survey (iFurusawa et al.1 12008 ) and JK imagin g from 
the UKIRT IR Deep Sky Survey (UKIDSS, .Lawrence et all 
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20071): for thes e data, we utilized the TiT-selected catalog from 
Wilhams et alj (120091) . The cluster field also has deep Spitzer 
imaging available in the four IRAC bands (3.6-8.0/im) and 
MIPS 24/im as part of the Spitzer public legacy surveyof the 
UKIDSS Ultra Deep Survey (SpUDS, PI: J. DunlopQ We 
matche d the JT-band selecte d catalog to the IRAC data fol- 
lowing [^govichetal] (120061) . 

CIG J02 18.3-05 101 was originally detected as a strong 
overdensity of galaxies at 1.5 < Zphot < 1-7 with a mean 
surface density > 20a higher than that of galaxies in the 
same redshift range in the UKIDSS Ultra-Deep Survey 
(iLawrence et aL 2007). Follow-up spectroscopy confirmed 10 
members that define a narrow redshift peak at 1.622 < z < 
1.634 with additional mem bers at z = 1.609 and z = 1.649 
(PaplO: 'Ta naka et"al]l2010l) . In our analysis, we include 80 
more members that are s elected from ph otometric redshifts 
determined with EAZY (Brammeret al. 2008): we use the 
limit for the integrated redshift probability defined by the 
spectroscopically confirmed members (T'. > 0.3: see PaplO 
for a detailed description of the photometric selection). 

Analysis of XMM-Newton data in this field also reveals a 
weak detection consistent with extended emission from the 
cluster. The cluster mass estimated from its X-ray luminosity 
is ^ 1.1 X lO''* Mq and consistent with its estimated dynam- 
ical mass of ~ 4 X 10''* M0. To define the cluster center, we 
select the massive (spectroscopically confirmed) member lo- 
cated at the peak of the X-ray emission: its J2000 coordinates 
are (2:18:21.07, -5:10:32.84), and all of the cluster galaxies 
lie within Rpmj ^ 1 Mpc of this member. 

2.1. Spectral Energy Distributions 

We fit the 10-band galaxy photometry (0.4/im< Xot>s < 
8/im|3 with the 2007 version of the Bruzual & CharloJ 
(|2003) stellar population synthesis models u sing a Chabrier 
initial mass function (for more details, see iPapovich et alj 
120061) . We find that models with Solar metallicity best re- 
produce the rest-frame colors and scatter of the cluster red se- 
quence galaxies (PaplO), and we allow the models to range in 
age from lO*" to 2 x 10' ° yr: we also include dust attenuation 
using the Calzetti et alj (12000.) law with color excess values 
ranging from £(5-^^) = 0.0-0.7. We allow for a range of star 
formation histories parameterized as a decaying exponential 
with an e-folding time r, where at any age t the star forma- 
tion rate is 'i'(t) ^ exp(-f /r) and r ranges from 1 Myr to 100 
Gyr (corresponding approximately to instantaneous bursts to 
constant star formation, respectively). In the model fitting, 
we add a <j/fi, = 5% error in quadrature to the photometric 
uncertainties on each band to account for mismatches in the 
multi-band photometry, and for the fact that the models do not 
continuously sample the model parameter space. Our sample 
of 92 cluster members has an average best-fit age of ^ 500 
Myr and an average best-fit r of ^ 500 Myr: 35 members 
have SED-derived star formation rates >5 Mq yr"'. 

2.2. MIPS 24fim Fluxes 

The cluster field was imaged with Spitzer MIPS as part of 
the legacy UKIDSS Ultra Deep Survey (SpUDS PI: J. Dun- 
lop). We extracted sources from t he public SpUDS 24/im map 
using StarFinder (iDiolaiti et al.1 IIOOO.) . an IDL-based PSF- 
fitting code designed for crowded fields. The 24/tm image 
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Fig. 1. — Rest-frame (U -V)ab color vs. stellar mass deteiTnined from fit- 
ting spectral energy distributions (SEDs; assuming Chabrier IMF) to the 10- 
band photometry: shown are the measured (top) and dust-corrected (bottom) 
color-mass diagrams. Open circles d enote members selected using photo- 
metric redshifts from EAZY t'Brammer et al. 2008) and filled circles spectro- 
scopically confirmed members from Papovich et al. (2010) and Tanaka et ;3] 
(2QI0); the latter tend to be on the blue edge because the spectroscopy favors 
members with emission lines. The 17 cluster members detected at 24/^m 
are marked with open stars; note the number of IR luminous members that 
remain red even after being corrected for dust extinction. The CIG J0218.3- 
0510 members have a color distribution similar to that observed in the field 
at z ~ 2 I Brammer et al. 2009 ) and, in contrast to galaxy clusters at z < 1 , the 
members populate a blue sequence and span a range in color. 

has 1.245 "pixel"', and we derived a model PSF and aper- 
ture corrections from the brightest isolated sources from the 
SpUDS image. The catalog includes all sources detected with 
S/N > 5; this corresponds to a flux of ^ 40/iJy. Using simu- 
lated sources based on the PSF and injected into the map, we 
determine that the catalog is 80% complete at this flux level. 
The measured 24/im flux es are converted into total infrared 
luminosities (L/r) using the lCharv & Elba j (12001) templates. 
Recent Herschel studies indicate that while this technique is 
very accurate at z < 1.5, extrapolations from monochromatic 
24/im fluxes overestimate the true L/« by factors of 2 - 7 at 
z > 1.5 (e.g. Nordon et al. 2010). Finally, star form ation rates 
are calculated from L/r using the prescription of iKennicuttI 
(IT998) . adjusted to the Chabrier IMF 

3. STELLAR MASSES, STAR FORMATION RATES, 
AND ENVIRONMENT 



By fitting SEDs to the 10-band photometry (0.4/tm< Ao/„ < 
8/im), we are able to measure accurate rest-frame colors (AB 
" http://ssc.spitzercaltech.edu/spitzermission/observingprograms/legacy/spuds/ system) and Stellar masses as well as correct for dust extinc- 
tion. Figure [T] shows the measured and dust-corrected rest- 
frame {U - V) color versus stellar mass for the 92 cluster 



At z = 1.62, this con'esponds to a rest-frame wavelength range of 
0.15^m< Xresi < 3.0/im. 



Reversal of the Star Formation-Density Relation in a Cluster at z = 1 .62 
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Fig. 2. — Star formation rate vs. stellar mass for 24/im detected cluster 
galaxies at z = 1 .62; symbols are as in Fig.[T] The most strongly star-forming 
systems are also some of the most massive cluster members. These galaxies 
follow the same trend that is observed in field galaxies at z ~ 2 (fitted relation 
and upper semi-interquai'tile range of 0.16 dex shown as solid and dotted 
lines, respectively: iDaddi et alj20Q7i) . All 17 of the 24/im detected members 
have hm > 10" Lq, in stark contrast to the relative dearth of such massive, 
IR-bright galaxies in clusters at z < 1 . 

galaxies within Rgro / ^ 1 M pc of the BCG. As demonstrated 
by, e.g. IWvder et al.l (120071) . the galaxies begin to separate 
into the well-known bimodal distribution only when the col- 
ors are corrected for extinction. 

The CIG J0218.3-0510 members dif fer from their coun- 
terparts in cluster s at z < 1.2 (e.g. iHolden et alJ 120061 : 
iRettura et al.ll2010l) in that they have a color-stellar mass dis- 
tribution that is surprisingly simila r to that observed in the 
field at z - 2 dBrammer et alJ '20Q9): the CIG J0218.3-0510 
members define a strong blue sequence and span a range 
in (dust-corrected) color (Fig. 2, bottom panel). While 
these members populate a red sequence in the (z-J) color- 
magnitude diagram (Pap 10), the correlation between rest- 
frame iU-V) color and stellar mass is visibly weakeQ How- 
ever, the most massive galaxies (M, ^ 10" Mq) are still the 
reddest, i.e. they have the oldest stellar populations. 

Seventeen of the members are detected at 24/im (Fig. [U, 
and the most IR luminous members include some of the most 
massive cluster galaxies (Fig.|2]i. Figure |2] shows that these 
IR luminous members follow the same trend of increasing 
star form ation with stellar mass that is observed in the field at 
z ~ 2 bvlPadd i et alj ( 12007 ,). Note that because we can only 
detect members that are IR-bright (L/s > 3 x 10" Lq), we are 
sensitive only to the upper envelope of the trend observed at 
z-2. 

Three of the CIG J0218.3-0510 members are Ultra- 
Luminous Infrared Galaxies (ULIRGs; L/s > 10'^ L©; see 
Fig. IS and the remaining 14 are LIRGs (10" <L/r < 
10'^ L0). In stark contrast, of the > 2000 galaxies in clusters 
at z < 1 studied with w ide-field (R„ro ; ^, 1 Mpc) mid-I R imag- 
ing (e .g. lGeach et al.i 2006; Sa intonge et al. 2008; Smit h et all 
120 101) . only one is a ULIRG and it lies outside the core of 
the Bullet Cluster, a well-known cluster-cluster merger at 

'^ The astute reader will notice that the red sequence is tighter in (z — J) 
color because these filters correspond approximately to rest-frame (U —B) at 
z=1.62. 
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Fig. 3. — Relative fraction of stai'-forming cluster galaxies vs. local galaxy 
density for 24fiva detected members with stai' fomiation rates > 40 Mq yr"' 
(top; open stars) and members with SED-derived rates >5 Mq yr"' (bottom; 
open triangles). The open circles in both panels show the 24/^m detected 
members in the massive galaxy cluster CL 1358+62 at z = 0.33 (Tran et al] 
2009) and represents the well-established trend documented at z < 1 ; these 
circles ai'e offset slightly in log E for claiity. Errorbai's are determined from 
bootstrapping the data in each bin 1000 times. The fraction of star-forming 
members is highest in the regions with the highest galaxy density, a reversal 
of the trend at z< \. Alz= 1.62, we have reached the epoch when massive 
cluster galaxies are still forming a significant number of new stai's. 



z = 0.297 ( Chung et al.ll20Tol) . The higher fraction of IR lu- 
minous galaxies at z = 1 .62 is likely to be driven by the over- 
all evolution o f the IR luminosity function in clusters (e.g. 
iBai et al.ll2009l) : however, more deep IR imaging of clusters 
at z > 1 is needed to confirm this trend. 

In our analysis, we assume that the 24fim sources are dom- 
inated by star formation and do not harbor active galactic nu- 
clei (AGN). From the SED fits, we find that two of the 24/im 
detected members do have emission at 8.0/xm (IRAC chan- 
nel 4) that deviates from the stellar fit, i.e. have a power-law 
component indicative of an AGN, and one of these galaxies 
is the most IR luminous member with log(L/R) [L0]~ 12.3 
(see Fig. |2]l. However, studies find that most of the emis- 
sion( > 70%) in ULIRGs is due to star formation (Farrah et al] 
I2OO8I) . thus these members are likely to have both strong star 
formation and an AGN component. Note that at z = 1 .62, the 
7.7/im PAH band lies partly in the 24/im channel. Because 
part of the IR luminosity is due to star formation, we include 
both 24/im members in our analysis; repeating our analysis 
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without these two members confirms that our overall results 
do not change. 

In the cluster core (Rproj = 0.5 Mpc), the star forma- 
tion rate density from the 24/im photometry alone is ^ 
1700 Mq yr"' Mpc"^; we stress that this is likely a lower limit 
given the 24/im imaging cannot detect any members with 
SFR/R < 40 M0 yr"', i.e. with star formation rates typical 
for IR-detected galaxies in clusters at z < 1 . Only one other 
galaxy cluster at z = 1.46 has a comparably high star forma - 
tion rate in its core C Havashi et al.| !2010: Hi lton et"al]|2010l) . 
In comparison, studies of IR-detected galaxies in clusters at 
z < 1 find that star-f orming members are strongly segregated 
a.tRf,.„i > 0.5 Mpc (iGeach et alJ l2006t ISaintonge et"aLll2008l: 
iKovama et alll2008h . 

Given the high star formation rate in its core, does CIG 
J0218. 3-0510follow the well-established trend atz< 1 ofde- 



creasing star-for mation with increasing galaxy density (e.g. 
Hashimoto et"an [1991 lEUingson et all l200!t iGomez et all 



20031) ? FigurelSlcompares the relative fraction of star-forming 



members to passive members as a function of local galaxy 
density ( S) which is defi ned by distance to the 10th nearest 
neighbor (Dresslei^'1980). We use star-formation rates de- 
rived from the 24/im imaging as well as from the SED fit- 
ting because the two independent star formation tracers com- 
plement each other and provide an important check of our 
results: The 24/im imaging is a robust measure of the dust- 
obscured star formation but only detects the most active mem- 
bers (> 40 M0 yr"') while the SED fitting is sensitive to lower 
levels of unobscured activity (>5 Mq yr"'). 

Both star formation tracers confirm that the relative frac- 
tion of active members is highest in the regions of highest 
galaxy density (Fig. O, i.e. exactly opposite to that ob- 
served in clusters at z < 1. A Spearman rank test supports 
with > 97% confidence (> 2cr significance) that the relative 
fraction of IR luminous members increases with increasing 
galaxy density from ~ 8% at E -- 10 gal Mpc"^ to -- 25% at 
E > 100 gal Mpc"^. We stress that excluding the two candi- 
date AGN does not change the trend, and the robustness of 
this result is underscored by the fact that we see the same 
trend using the SED-derived rates. While studies of galax- 
ies in the field at z ~ 1 (E < 10 gal Mpc"^) find that the 
star for mation-de nsity relation is beginning to turn over at this 
epoch (lElbaz et al. 2007; Cooper et al. 2008) . this is the first 
confirmation of such a reversal in the significantly higher den- 
sity environment of galaxy clusters. 

The measured IR luminosities correspond to specific star 
formation rates (star formation rate divided by stellar mass; 
SSFR) per Gyr of ^ 1 -20: these active members can more 
than double their stellar masses in the next Gyr (by z ^ 
1.2). However, to reproduce the relatively homogeneous stel- 
lar ages measured in massive cluster galaxies at z < 1 (e.g. 



iBlakesIee et al.ll2006l:PRan"etani2007HMei et al.ll2009h . these 
IR luminous members cannot maintain such a high SSFR for 
even a Gyr. The current star formation must be quenched 



rapidly and any later bursts of activity cannot add a substan- 
tial amount of new stars, at least not in the massive members 
(log(M,) [M0]> 10.6) that must populate a well-defined red 
sequence by z '^ 0.8. 

4. CONCLUSIONS 

We measure the rest-frame colors (dust-corrected), IR lu- 
minosities, star formation rates, and stellar masses of galaxies 
in CIG J0218. 3-0510, a Spitzer-selected cluster at z = 1.62, 
by fitting spectral energy distributions to photometry in 10 
bands (0.4/im< \„bs < 8/i m) and with deep 24/im imaging. 
The cluster sample (PaplO; lTanakaet"al]l201 0) is composed 
of 12 spectroscopically confirmed members and 80 members 
selected from photometric redshifts measured using EAZY 
(Brammer et al. 2008); all members are within Rpmj ^ 1 Mpc 
of the massive cluster galaxy located at the peak of the X-ray 
emission. 

The 92 cluster members have a color-stellar mass distribu- 
tion that is surprisingly similar to that observed in field galax- 
ies at z ~ 2. When corrected for dust, the cluster members 
define a strong blue sequence and span a range in color, indi- 
cating a substantial amount of recent and ongoing star forma- 
tion in the cluster core. This dramatic level of activity is un- 
derscored by the 17 members detected at 24/im. In the cluster 
core (Rproj < 0.5 Mpc), the star formation rate density from 
the IR luminous members alone is ^ 1700 Mq yr"' Mpc"^; 
the true value is likely to be higher given that we only in- 
clude members with SFR/r > 40 Mq yr"' . These IR luminous 
members also follow the same trend of increasing star forma- 
tion with stellar mass that is observed in the field at z ^ 2. 

We discover the striking result that the the relative frac- 
tion of star-forming galaxies increases with increasing local 
galaxy density in CIG J0218. 3-0510, a reversal of the well- 
established trend at lower redshifts and in line with recent 
work at z ^ 1 .46 that suggests enhanced star formation in 
cluster cores. Measurements using star formation rates de- 
rived from the 24/im imaging and from the SED fitting pro- 
vide independent confirmation that the relative fraction of 
star-forming galaxies triples from the lowest to highest den- 
sity regions. By pushing into the redshift desert (z > 1 .6), we 
are able to reach the epoch when massive cluster galaxies are 
still forming a significant number of their stars. 
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